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ABSTRACT 

We examine the chemical properties of five cosmological hydrodynamical simulations 
of an M33-like disc galaxy which have been shown previously to be consistent with 
the morphological characteristics and bulk scaling relations expected of late-type spi- 
rals. These simulations are part of the Making Galaxies In a Cosmological Context 
(MaGICC) Project, in which stellar feedback is tuned to match the stellar mass - halo 
mass relationship. Each realisation employed identical initial conditions and assembly 
histories, but differed from one another in their underlying baryonic physics prescrip- 
tions, including (a) the efficiency with which each supernova energy couples to the 
surrounding interstellar medium, (b) the impact of feedback associated with massive 
star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling 
of Type II supernovae remnants, (d) the treatment of metal diffusion, and (e) varying 
the initial mass function. Our analysis focusses on the resulting stellar metallicity dis- 
tribution functions (MDFs) in each simulated (analogous) 'solar neighbourhood' (2—3 
disc scalelengths from the galactic centre) and central 'bulge' region. We compare 
and contrast the simulated MDFs' skewness, kurtosis, and dispersion (inter-quartile, 
inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empiri- 
cal solar neighbourhood MDF and local group dwarf galxies. We find that the MDFs 
of the simulated discs are more negatively skewed, with higher kurtosis, than those 
observed locally in the Milky Way and local group dwarfs. We can trace this difference 
to the simulations' very tight and correlated age-metallicity relations (compared with 
that of the Milky Way's solar neighbourhood), suggesting that these relations within 
'dwarf discs might be steeper than in L* discs (consistent with the simulations' star 
formation histories and extant empirical data) and/or the degree of stellar orbital re- 
distribution and migration inferred locally has not been captured in their entirety, at 
the resolution of our simulations. The important role of metal diffusion in ameliorating 
the over-production of extremely metal-poor stars is highlighted. 
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1 INTRODUCTION 

The relative number of stars of a given metallicity in a given 
environment, whether it be the local stellar disc, central 



spheroid/bulge, and or baryonic halo - the so-called metal- 
licity distribution function (MDF) - has embedded within 
it, the time evolution of a system's star formation, assem- 
bly/infall, and outflow history, all convolved with the ini- 
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tial mass function (IMF) |Tinslevlll980l ). Seminal reviews of 
the di agno stic powe r of th e MDF include those of iHavwoodl 
l|200ll ) and lCaimmil l|200Sl ). 

Well in advance of our now empirical appreciation of 
(a) the hierarchical assembly of galaxies from sub-galactic 
units, (b) the ongoing infall of fres h material from hal os 
to discs (e.g. High- Velocity Clouds: I Gibson et al.l l|200ll )). 
and (c) the ongoing outflow of enriched material from discs 
via stellar- and supernovae -driven winds/fountains (e.g. 
iMcClure-Griffiths et al.ll2006l l . it was recognised that the lo- 
cal MDF provided crucial evidence that the Milky Way (and 
presumably galaxies as a whole ) did not behave as a 'cl osed- 
box', in an evolutionary sense (jPagel fc Patchettlll975l ). 

This latter recognition was perhaps best m anifest in 
what became known as the 'G-dwarf Problem' l|Hartwickl 
1 19761 ); specifically, a simple model in which gas was not al- 
lowed to infall or outflow from the system would necessarily 
lead to a significant population of long-lived, low metallic- 
ity, stars in the solar neighbourhood, with ~20% of the stars 
locally predict ed to possess metallicities below [Fe/H]«— 1 
l|Tinslevl [l980). In nature, such a population is not ob- 
served, with the empirical fraction of local low-metallicity 
stars being ~2 orders of magnitude s maller than the afore- 
mentioned closed-box p redictions fe.g. lKotoneva et al ] |2002l ; 
ICasagrande et al1l201lh . 

Since this recognition of its fundamental importance, 
the MDF has acted as one of the primary constraints 
/ boundary conditions against whic h all analytical (e.g. 
Schorck et al]|2009l; iKirbv et alj|201ll). semi-n umerical (e.g. 
Chiappini et al . 2001; Fcnncr & Gibson 2003), and chemo- 



dynamical (e.g. iRoskar et alj 2008: Sanchcz-Blaz quez et alj 
120091 : iTissera et alj l201ll : ICalura et alj l2012h models are 
compared. 

From a chemo-dynamical perspective, recent work has 
focused on the sensitivity of global metal re-distribution to 
different physi cal prescriptions , with in the context of the 
OWLS project Iwiersma et alj (|201ll ); at higher redshift, a 
similar, equally comprehensive , study was undertaken by 
ISommer-Larsen fc Fvnbol (|2008l ). In both cases, the empha- 
sis was placed on the whereabouts of the 'missing metals' - 
i.e., metals thought to reside in the Warm-Hot Intergalac- 
tic Medium (WHIM) and/or halos of massive galaxies, but 
have thus far proven challenging to detect directly^ While 
not fully cosmological, the reader is also referred to the 
excell ent chemo-dynamical work of iKobavashi fc Nakasatol 
(2011), for a complementary analysis of a simulated Milky 
Way-like system. 

Each of the above chemo-dynamical studies examine 
cursory aspects of the MDF 'constraint', but the focus for 
each was never meant to be a comprehensive analysis of the 
dispersion and higher-order moments of the shape charac- 
teristicsQ nor their link to the associated age-metallicity re- 
lations, star formation histories, and putative G-dwarf prob- 
lem; such higher-order moments include the MDF skewness, 



kurtosis, and inter-quartile, inter-decile, inter-centile, and 
inter-tenth-percentile regions. 

The skewness of an MDF can be a reflection of both 
the classical G-dwarf problem and the slope of the age- 
metallicity relation (AMR); kurtosis is often thought of as 
being a measure of the 'peakedness' of the MDF (e.g., by 
how much the peak is 'flatter' or 'peakier' than a Gaus- 
sian), while in practice it is often more sensitive to the pres- 
ence of 'heavy' tails, rather than the shape of the peak; the 
inter-quartile, -decile, etc., regions probe both the effects 
of star formation histories and AMRs and, in the case of 
the inter-centile and (especially) the inter-tenth-percentile 
regions, the impact of metal diffusion on the extreme metal- 
poor tail of the distribution. In the context of cosmological 
chemo-dynamical disc simulations, to our knowledge, our's 
is the first quantitative discussion of these higher-order mo- 
ments of the MDF. 

Further, from an observational perspective, the re- 
cent re-calib rations of the orig i nal Ge neva- Copenhagen Sur- 



1 iTumlinson et all 1120111 ) is an excellent example of recent efforts, 
though, to characterise the properties of these difficult-to-observe 
bary on reservoirs . 

2 cf. IKirbv et al. I J201lll . though, for a study of the higher-order 
moments of the MDFs of Local Group dwarf spheroidals which is 
similar in spirit to our work here on disc galaxies. 



yey ( GCS: iNordstrom et alj d2004ht by iHolmberg et all 
l|2009l ) and ICasagrande et alj (|201 lh has made for a 
timely investigation of the predicted characteristics of the 
MDFs of simulated disc galaxies. Parallel developments 
slightly further afielojf] include targeted MDF studies of the 
thin— thick disc transitio n region and the thi ck disk proper 
JSchlesinger et alj I2012T ), the stellar halo (pc horck et al 
| 2009h . and the Galactic bulge (jBensbv et"aill201ll : lHill et al 

This paper aims to fill a gap in the literature, by mak- 
ing use of a new suite of fully cosmological chemo-dynamical 
simulations whose properties have been shown to be in re- 
markable agreement with the basic scal ing laws to which 
late-type disc gala xies adhere in nature (|Brook et al J 12013 ; 
iMaccio et alj|20l3 ). The simulations themselves are outlined 
briefly in [J2] alongside a description of the adopted anal- 
ogous 'solar neighbourhood' regions. The associated age- 
metallicity relations (AMRs) are presented in SJ3] the need 
for this will become apparent when analysing the higher- 
order moments of the MDFs within these regions and, in 
particular, their metal-poor tails (Q. Our results will then 
be summarised in 35] 



2 SIMULATIONS 

In what follows, we analyse five cosmological zoom vari- 
ants of the 'scaled-do wn' M33-like disc galaxy simulation 
(g 15784) described by IfBrook et alj (12013 ). The initial con- 
ditions are identical fo r each realisati o n, an d taken from the 
eponymous R15784 oflStinson et alj (|2010h after re-scaling 
(e.g. IKannan et alj I20l3 ) the mass (length) scales by a 
factor of eight (two). Differences in the underlying power 
spectrum that result from this re-scaling are minor (e.g . 
Springel et al.ll200Sl;lMacci6 et alj|20l3 ; IKannan et aLll20l3 ; 
Vera-Ciro et al.l 12013 )7 and do not affect our results. The 
virial mass of the scaled gl5784 is 2xlO n M Q , with ~10 7 
particles within the virial radius at 2=0, with a mean stel- 
lar particle mass of ~6400 Mq . A gravitational softening of 



3 Spatially speaking, in relation to that of the solar neighbour- 
hood region probed by the GCS. 
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e=155 pc was used; to ensure that gas resolves the Jeans 
mass, rather than undergoin g artifical fragmentation, pres - 
sure is added to the gas, after lRobertson fc Kravtsovl (|2008f l . 
Further, a maximum density limit is imposed by setting the 
minimum SPH smoothing length to be 1/4 that of the soft- 
ening length 

Each of the five simulations was evolved using the 
gravitational N-body + sm oothed particle hydr odynam- 
ics (SPH) code Gasoline (|Wadslev et all 12004 ). Metal- 
dependent cooling of the gas, unde r the assumption of ion- 
isation equilibrum, is applied , after IShen et all (I2010). cou- 
pled to a uniform, evolving, lHaardt fc Madaul |l99rj ) ion- 
ising ultraviolet background. Our re ference/fiduc i al sim ula- 
tion ( 1 lmKroupa) was introduced bv lBrook et al.l (|2012f ). in 
the context of its outflow and angular momentum character- 
istics. The structural and kinematic properties (e.g., rotation 
curves, bulge-to-disc decomposition, ratio of rotational-to- 
anisotropic support, etc.) of the simulations pr esented here 
are in distinguishable from those presented in iBrook et al.l 
(2012), to which the reader is referred for supplementary 
details. 

When gas reaches a sufficiently cool temperature - 
T< 10,000— 15,000 K - and resides within a sufficiently dense 
environment - n t h>9.3 cm 
stars according to ^f- =c *-jf^ 
tion efficiencylf] At is the timestep between star formation 
events (0.8 Myrs, here), M ga s is the SPH particle mass, td yn 
is the SPH particle's dynamical time, and AM* is the mass 
of the star particle formed. 

We have extended the chemical 'network' of Gaso- 
line from oxygen and iron, to now also track the evo- 
lution of c arbon, nitrog e n, ne on, magnesium, and sili- 
con. After Raiteri et al.l l| 19961 ). power law fits to the 



! -0 it becomes eligible to form 
where c* is the star forma- 



IWooslev fc Weaver! i|l995h Z=0.02 SNell yields were gen- 
erated for the dominant isotopes for each of these seven ele- 
ments; a further extension was implemented, in order to in- 
clude the Ivan den Hoek fc Groenewegenl (|l997l ) metallicity- 
dependent carbon, nitrogen, and oxygen yields from asymp- 
totic giant branch (AGB) stars. By expanding upon the 
chemical species being tracked, the earlier concern regard- 
ing the underprediction of the global metallicity by a fac- 
tor of ^2 (and the consequent underestimate to the SPH 
cooling and star form ation rates) is naturally alleviated 
l|Pilkington et al.ll2012T ). We note in passing that all abun- 
dances (and rat ios) presented here ar e relative to the solar 
scale defined bv lAsplund et al.l (|2009h . 

Feedbac k from supern o vae (S Ne) follows the blastwave 
formalism of lStinson et al.1 l|2006h . with 100% of the energy 
(10 51 erg/SN) thermally coupled to the surrounding ISM. 
Cooling is disabled for particles within the blast region (cor- 
responding to the radius of the remnant when the interior 



4 In comparison, the original IStinson et al.l (|2010h simulations 
used a minimum SPH smoothing length of e/100, resulting a dra- 
matic increase in computational time, but with only minimal im- 
pact on the simulation itself. 

5 The star formation density threshold n t ^ corresponds to 
the maximum density gas can reach using gravity — i.e., 
n th =32M gas /e 3 . 

6 The star formation efficiency c* was taken to be 10% for all the 
runs, except for llmChab, for a value of 7.5% was adopted. 



pressure has been reduced to that of the pressure of the am- 
bient ISM) for a timescale corresponding to that required 
to cool the hot interior gas to T~10 4 K(j Bearing in mind 
the 0.8 Myr timesteps of our runs, we impose a minimum 
cooling 'shut-off time' which matches this value0 

We employ the "MaGICC" (Making Galaxies In 
a Cosmolog i cal C ont ext) feedback mode l described by 
IBrook et alT(|2012l i and lStinson et al.1 (|2012l i. taking into ac- 
count the effe ct of energy feedback from massive stars into 
the ISM3 (cf. iHopkins et ai] (|201lT )1. While a typical mas- 
sive star might emit ~10 53 erg of radiation energy during 
its pre-SN lifetime, these photons do not couple efficiently 
to the surrounding ISM; as such, we only inject 10% of this 
energy in the form of thermal energy into the surrounding 
gas, and cooling is not disabled for this form of energy input. 
Of this injected energy typically 90-100% is radiated away 
within a single dynamical time. 

The defa u lt ini tial mass function (IMF) is that of 
iKroupa et all (|l993h ; the llmChab run incorpo rates the 
more contemporary (and currently favoured) IChabrierl 
(2003) functional form; per stellar generation, the latter pos- 
sesses a factor of ~4x the number of SNell as that of the for- 
mer. Finally, the t reatment of metal diffusion within Gaso- 
line is detailed bv lShenet all |2010f ); a diffusion coefficient 
C=0.05 has been adopted for our runs, except for one simu- 
lation for which diffusion was prohibited (HmNoDif f )0 The 
primary numerical characteristics of the five simulations em- 
ployed here are listed in Table [T] 

For our MDF and AMR analyses, for each simulation we 
identify an analogous region to that of the Milky Way's 'solar 
neighbourhood', defined to be a radial range from 3.0 to 3.5 
disc scalelengths (see Table[TJ) and to lie within 500 pc of the 
galactic mid-plane. The fraction of accreted stars in these 
high-feedback runs is negligble; as such their contamination 
in the 'solar neighbourhood' is equally negligible. Conse- 
quently, there was no need to undertake the sort of kinematic 
decomposition of the orbital circularity ej = J z / Jdrc(E) 
distributiorl"! that was needed to isolate disc/in-situ stars 
from spheroid/accreted stars in our parallel analysis of the 
MDFs of the more massive (and accretion-contam inated) 
IStinson et al.l (|2010h simulations l|Calura et al.l l2012)FI 



To use the terminology of iGibsonl lll994) , the relevant radius 
and timescale correspond to R mer gc and t coo i, respectively. 

8 Save, for the one run for which this restriction was relaxed 
(HmNoMinShut). 

9 Except for the one run included here without radiation energy 
(HmNoRad). 

10 Our 'rin rl i+Ti icir^n ' run nnccocooa A/TT^Th 1 nnrl n\i «tti i r 1 n 1 1 r-\~\ a r a ft rn- 



llGovernato et al.l 


20101). the latter 


ysis was shown in 


Pilkington et al.l 



(2012). This similarity can be traced to the less efficicient metal 
diffusion adopted for the DG1 runs (i.e., C =0.01 vs th e C=0 .05 
now employed for our Gasoline runs, after IShen et al] (l2010h i. 

11 Where J z is the z-component of the specific angular momen- 
tum and J c irc(E) is the angular momentum of a circular orbit at 
a given specific binding energy. 

12 Note, this was confirmed by undertaking a kinematic decom- 
position of 1 l mKro upa using the modi fied technique introd uced by 
lAbadi et al.l d2003f) . and employed bv lCalura et al. (2012); specif- 
ically, none of our conclusions were contingent upon the need for 
a kinematic decomposition. More quantitatively, only ~3% of the 
stars in our simulated 'solar neighbourhoods' would be kinemat- 
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Galaxy 


IMF 


c* 


eSN 


SR 




Stellar Mass 


Scale Length 


Vertical Gradient 


Radial Gradient 


llmKroupa 


Kroupa 


0.1 


100% 


10% 


15000 


7.1xl0 9 


2.34 


-0.064 


-0.012 


llmChab 


Chabrier 


0.075 


100% 


10% 


10000 


1.3xl0 9 


2.78 


-0.017 


-0.026 


llmNoRad 


Kroupa 


0.1 


100% 


0% 


15000 


9.1xl0 9 


1.58 


-0.027 


-0.045 


HmNoMinShut 


Kroupa 


0.1 


100% 


10% 


15000 


14.0x10 s 


1.71 


-0.008 


-0.020 


HmNoDiff 


Kroupa 


0.1 


100% 


10% 


10000 


2.1xl0 9 


1.43 


-0.013 


-0.028 



Table 1. Primary para meters employed for the five sim ulations analysed in this work. Column (1): simulation/galaxy name; Column (2): 
adopted IMF (Kroupa ^Kroupa et al.l l|l993h ; Chabrier ^Chabrier] l|2003h : Column (3): star formation efficiency; Column (4): thermalised 
SN energy fraction coupled to the ISM; Column (5): thermalised massive star radiation energy fraction coupled to the ISM; Column (6): 
maximum allowable gas temperature for star formation; Column (7): present-day stellar mass (in solar masses) within the virial radius; 
Column (8): stellar disc exponential scalelength (in kpc); Column (9): vertical [Fc/H] gradient (in dcx/kpc); Column (10): radial [Fe/H] 
gradient (in dex/kpc). 



We first show the inferred star formation histories 
(SFHs) of the solar neighbourhoods associated with each of 
the five simulations (Fig[l|. Several important points should 
be made, before analysing the AMRs and MDFs. Qualita- 
tively speaking, the SFHs of these regions within llmKroupa, 
llmNoMinShut, and llmNoRad are similar to those seen in 
gas-rich dwarfs like NGC 6822, Sextans A, W LM, and 
to some extent, the LMC jDolphin et all 120051 ). In that 
sense, they are (not surprisingly) different from the typi- 
cal exponentially-decaying SFH (timescales of ~5— 7 Gyrs) 
inferred for the Milky Way's solar neighbourhood (e.g. 
iRenda et al"1l2005l ), and so we should not expect identical 
trends in the ancillary AMRs and MDFs, as those observed 
locally. Indeed, we will show this to be case momentarily, 
but our interest here is more in identifying trends, rather 
than exact star-by-star comparisons. 

The one simulation which shows an exponentially- 
declining SFH at later times is that of HmNoDiff ; the lack 
of diffusion here acts to minimise the 'spread' of metals to 
a degree that star formation is restricted (preferentially) to 
much less enriched SPH particles (in part, because the cool- 
ing then becomes less efficient for a greater number of SPH 
particles, which has a greater impact at later times where 
there are fewer efficiently cooling metal-enriched SPH par- 
ticles out of which to potentially form stars. We will return 
to the special case of the 'no diffusion' model shortly. 

The SFH of llmChab also shows a distinct behaviour 
relative to the llmKroupa fiducial. Specifically, it is sig- 
nificantly lower, and relatively constant, at all times; in 
spirit, this is similar to the inferred SFH of the LMC (e.g. 
[Holt zman et al.lll999l ). This is reflected in the stellar mass 
at z—0 being significantly lower than llmKroupa, which in 
turn aids considerably in bringing its properties into close 
agreement with es sentially all traditional scaling relations 
iBrook et al.1 (|2012l ). This behaviour is driven by (a) the fac- 
tor of four increase in the SNe per stellar generation (via the 



ically classified as 'bulge/spheroid' stars, impacting on the vari- 
ous MDF metrics to be discussed later at the <3% level (smaller 
than the uncertainty associated with the treatment of extreme 
(>5<t) outliers - see §4). In light of this negligible impact, we 
have avoided imposing any personal preferred kinematic decom- 
position scheme into the analysis. 




Time [Gyr] 

Figure 1. Star formation histories of the solar neighbourhoods 
associated with the five simulations; colour-coding is as noted in 
the inset to the panel. 

more massive star-biased IMF), and (b) the reduced maxi- 
mum temperature for star formation (as noted earlier). 

The subtle effect of allowing the minimum shut-off time 
for radiative cooling of SN remnants to become prohibitively 
small in high-density regions (in practice what this means is 
that the shut-off time becomes smaller than the timestep of 
0.8 Myrs) can be seen in the llmNoMinShut curve of FigQ] 
Specifically, SPH particles affected by this effectively cool 
'instantly' within the same timestep, without any delay. 
Hence, the particles in question become 'available' for star 
formation much sooner than they might otherwise; this has 
the effect of 'boosting' the star formation relative to that of 
the fiducial llmKroupa. 



3 AGE-METALLICITY RELATIONS 

As noted earlier, the MDF bears the imprint of a re- 
gion's star formation history (SFH), convolved with its age- 
metallicity relation (AMR). Having introduced the 'solar 
neighbourhood' SFHs in ij2] we now present their associated 
AMRs in Fig[2] The time evolution of the [Fe/H] abundances 
are shown for each of the five simulations listed in Table [1] 
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Colour-coding within each panel corresponds to stellar age, 
ranging from old (black/blue) to young (red). 

To provide a representative empirical dataset against 
which to compare, we make use of the recent re-calibration 
of the Geneva-Copen hagen Survey (CCS) presented bt 
Irlolmberg et alj (|2009l ) . The base CCS provides invaluable 
spectral parameters for ~ 17,000 F- and G- stars in the solar 
neighbourhood. Following iHolmberg et al.l (|2007l ). we define 
a 'cleaned' sub-sample by eliminating (i) binary stars, (ii) 
stars for which the uncertainty in age is >25%, (iii) stars 
for which the uncertainty in trigonometric parallax is >13%, 
and (iv) stars for which a 'null' entry was provided for any of 
the parallax, age, metallicity, or their associated uncertain- 
ties. The AMR for this 'cleaned' sub-sampled of ~4,000 stars 
is shown in the lower-right panel of Fig [2] A fifth criterion is 
applied for the determination of the highe r-order moments 
of the MDF shape; specifically, following IHolmberg et al.l 
|2007h and constructing an unbiased volume-limited sub- 
sample from the stars lying within 40 pc of the Sun. Doing 
so yields a smaller sample of only ~500 stars. While this 
does not impact on the shape characteristics of i|4] or the 
behaviour of the AMR, for clarity, we show the AMR in- 
ferred from the aforementioned sub-sample of ~4,000 stars 
in Fig [20 

It is worth re-e mphasising that we are using the 
IHolmberg et al.l (|2009h variant of the GCS solely as a useful 
'comparator' against which to contrast our various MDF 
metrics / higher-order moments. It should not be inter- 
preted as an endorsement of one solar neighbourhood MDF 
over another; there is a rich literature describing the various 
pros and cons of any number of potential selection biases 
within this (or any other) re-calibration of the G CS (e.g. 
ISchoenrich fc BinnevlliiOcI ; ICasagrande et alj|201ll ) and we 
are not equipped to enter into that particular debate. The 
GCS remains the standard-bearer for MDF analysis, reflect- 
ing the nature of (fairly) volume-limited and (fairly) com- 
plete nature, making it ideal for probing the active star form- 
ing component of the thin disc; other exquisite MDFs, in- 
cludi ng those of the aforem entioned ( predominantly) thick 
disc l|Schlesinger et alj|2012h and halo (jSchorck et al.ll2009l ) 
studies, are more suited for simulations targeting regions 
further from the mid-plane than we are doing here. Ide- 
ally, of course, we would like to replace the solar neigh- 
bourhood 'comparator' used here (the GCS) with an em- 
pirical sample more representative of star formation his- 
tories associated with massive dwarf spiral/irregulars (e.; 



ISkillman et"al]|2003l ; iDolphin et aLll2003l ; iKirbv et alj|201lj ) , 
but until the statistics, completeness, and accuracy of the 
age and metallicity determinations for such distance dwarfs 
reaches that of the solar neighbourhood, we are reluctant to 
compare (in detail) the predictions of the simulations with 
those of the observations. Having said that, we will comment 
on, in a qualitative sense, the AMR and MDF trends seen 
in our simulations and how they compare with said dwarfs. 

Several key points can be inferred from Fig [2] First, 
not surprisingly, the metallicities of the stars in the Milky 



13 The 'upturn' towards high-mctallicities at young ages in the 
GCS sample is likely traced to the very young Fm/Fp stars which 
are difficult to charac terise with Stromgren photometry alone 
llHolmbere et alj|2009h . 



Way's solar neighbourhood (GCS) are typically a factor of 
~5— 100 x higher at a given age compared with the five sim- 
ulations. This reflects the discussion of 321 in relation to the 
fact that the simulations in question are more similar to 
lower-luminosity disc galaxies (in terms of both mass and 
SFHs), rather than being Milky Way 'clones'. The simu- 
lations are consistent with the various s caling relations to 
which galaxies adhere (|Brook et al.ll2012f) ; as such, for their 
mass, their mean metallicities are a factor of ~3— 5x lower 
than that of the Milky WayE3 

More important for our purposes here, there are two ad- 
ditional characteristics which are readily apparent in FigO 
First, the AMR of the solar neighbourhood is essentially 
non-existent, save for a trace of old, metal-poor, stars. In 
contrast, the corresponding regions of the simulations show 
extremely correlated AMRs (especially those of the fidu- 
cial simulations, llmKroupa and llmChab). This is partly 
traced to the differences in the aforementioned SFHs, al- 
though the correlation persists (admittedly with larger scat- 
ter at a given age) even in HmNoDif f , the simulation whose 
SFH bears the closest resemblance to that of the Milky 
Way. The impact of these tightly-correlated AMRs mani- 
fest themselves significantly within the inferred MDFs, a 
point to which we will return in Qualitatively speak- 
ing, these tightly-correlated AMRs resemble those predicted 
by semi-numerical ga l actic chemical evolution models (e.g. 
Chiappini et al.ll200ll; iFenner fc Gibson! 120031 ; iRenda et al.l 



20051 ; iMolla fc DfadbOOsf T 



In the bottom right panel of Fig [2] we also overplot the 
AMRs inferred from the colour-magnitude diagram-derived 
star formation histor ies of the dwar f irregulars Sextans A 
l|Dolphin et all I2003T ) and IC 1613 l|Skillman et al.ll2003T l; 
like the Milky Way, neither are meant to be one-to-one 
matches to the 11m series of simulations, but in some sense 
they do provide a useful complementary constraint, in the 
sense that their respective star formation histories are not 
dissimilar to those shown in Fig [T] (in particular, those of 
llmKroupa, HmNoMinShut, and HmNoRad). Their associated 
AMRs, while lacking the statistics, completeness, and accu- 
racy of the GCS dataset necessary to make detailed quantia- 
tive comparisons, do show evidence of possessing somewhat 
stronger correlations. Again, the statistics of these dwarf 
systems' MDFs and AMRs make it difficult to say anything 
more regarding the degree of 'agreement' between the 11m 
series and that encountered in nature, but it is suggestive 
and certainly merits revisting once data comparable to that 
of the GCS becomes available for dwarf irregulars/spirals. 

Second, the scatter in [Fe/H] at a given stellar age is sig- 
nificantly smaller (compared with that of the Milky Way) 
in the three simulations where the injection of thermalised 
massive star radiation energy to the surrounding ISM is in- 
cluded (i.e., llmKroupa, llmChab, and HmNoMinShut). Ne- 



14 The MDFs and AMRs of systems more directly comparable to 
the Milky Way proper - i. e., the more massive 'parent' simulations 
to those employed here llStinson et al.l P2010l - are described by 
ICalura et alj d2012T> and Bailin et al. (2012, in prep), respectively. 
The significant contamination from accreted stars in these more 
massive simulations tends to impact upon both the scatter of the 
AMR and skewness/dispersion of the IMF, in a negative sense, 
relative to the high-feedback models here, for which the accreted 
fraction is negligible. 
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Figure 2. Age-metallicity relations (where metallicity=[Fe/H]) 
in the analogous solar neighbourhoods of the five simulations em- 
ployed here, in addition to the reference relationship found in 
the solar neigh bourhood of the Mil ky Way and t he dwarf irregu- 
lars S extans A jDolphin et al.ll2003l) and IC 1613 llSkillman et al.l 
2003). Colour-coding in each panel is by stellar age, ranging from 
black/blue (oldest) to red (youngest). 

glecting this feedback term, within the context of these cos- 
mological hydrodynamical disc simulations, acts to increase 
the scatter in [Fe/H], at a given in time, to a level compara- 
ble to that seen in Milky Way's solar neighbourhood!^ Not 
surprisingly, the one simulation for which metal diffusion 
was suppressed (HmNoDiff) possesses the largest scatter in 
[Fe/H] at a given age, particularly at early times/low metal- 
licities, where the neglect of diffusion is most problematic 
(again, a point to which we return in 



4 METALLICITY DISTRIBUTION 
FUNCTIONS 

Having been informed by the empirical and simulated so- 
lar neighbourhoods' SFHs and AMRs ($21 anci 33 we now 
present the [Fe/H] metallicity distribution functions (MDFs) 
for the same regions0 Fig [3] shows the MDFs (black his- 
tograms) for the five simulations, the Milky Way (GCS: 
lower right panel) and the Local Group dwarf Fornax (also, 



15 A secondary byproduct is also a mildly steeper radial abun- 
dance gradient, although the effect is minor - recall, Table [T] 

16 We confirmed that our conclusions are robust to the specific 
definition of the 'solar neighbourhood', by increasing its vertical 
range from ±0.5 kpc to ±2 kpc. Similarly, varying the radial range 
from 3.50±0.25 disc scalelengths, by ±1 scalelength has negligblc 
impact (recall from Table [T] that the metallicity gradients here 
are shallow). 



lower right panel, from iKirbv et al.l (|201ll )). The two sub- 
samples of the GCS are shown; in black, the aforementioned 
(33 sub-sample of ~4,000 stars (matching those shown in 
Fig El - i- e -i the 'cleaned' sub-sample, but without any dis- 
tance constraint applied, labeled 'GCS' in the lower-right 
panel), and in blue, the volume-limited sample (i.e., those 
lying within 40 pc of the Sun, labeled 'GCScut'). As stressed 
earlier, the shape characteristics of the GCS MDF are not 
contingent upon this latter cut; the labels 'GCS' and 'GCS- 
cut' will be employed to differentiate between the two, where 
relevant. Overlaid in each panel is simple 'best-fit' (single) 
Gaussian to the respective distributions (and their associ- 
ated full- width at half-maximum (FWHM) values). For the 
Forn ax dwarf, we use t he full sample of 675 stars taken 
from IKirbv et al.l (|201ll ). in order to show (perhaps) the 
best determined MDF for a representative local dwarf. Three 
caveats should be noted, in relation to the latter: (i) the sam- 
ple size is 1-2 orders of magnitude smaller than the GCS, 
not surprisingly, considering the challenging nature of this 
observational work; (ii) no analogous 'solar neighbourhood' 
can be identified within this dataset (it is simply all the stars 
in the sample covering a range of fields in Fornax); and (iii) 
the uncertainty in [Fe/H] for a given individual star in For- 
nax is ~0.5 dex, compared with the ~0.1 dex associated 
with individual stars in the GCS. Fornax is neither better 
nor worse than the GCS, as a comparator, so it is useful to 
at least show both, as they represent the state-of-the-art, 
observationally-speaking. 

Even before undertaking any quantitative analysis of 
the MDFs, it is readily apparent that the simulations (par- 
ticularly, llmKroupa, and llmChab) possess an excess of stars 
to the left (i.e., to the negative side) of the peak of the MDF, 
relative to the right, when compared with that of the GCS 
and Fornax (i.e., the simulated MDFs are more negatively 
skewed). This 'excess' of lower-metallicity stars are formed 
in situ during the first ~4 Gyrs of the simulations. The ex- 
ception to this trend is HmNoMinShut, for which the lack of 
significant star formation at early epochs (recall, Fig[T} and 
the extremely flat AMR at late times (Fig [2| conspires to 
present the narrow and symmetric MDF shown in Fig [3] As 
noted in iJU for both HmNoRad and HmNoDiff, the larger 
scatter in [Fe/H] at a given age manifests itself in the broader 
MDFs seen in Fig [3] 

It is worth delving deeper into the source of the broader 
MDF seen in, for example, HmNoRad, relative to the fidu- 
cial llmKroupa. Here, it is at high-redshift that the radiation 
energy has an impact on the regulation of star formation. 
HmNoRad has higher star formation at early times (Fig 1), 
but not at later times, primarily because it exhausts its avail- 
able gas, whereas with the radiation energy star formation 
is regulated during that crucial period when gas accretion is 
at its most active; this gas remains available at later times to 
form stars, resulting in the MDF of HmNoRad being broader 
relative to the fiducial. Ultimately, the length of time that 
gas spends in the disk before it forms stars shapes the MDF 
'width' here. With radiation energy included, this gas is in 
the disk for a longer period of time, meaning more metal 
mixing occurs. Linking back to the star formation histories 
of Fig 1, we note that most of the gas is accreted during 
the first ~6 Gyr, and one can see that the star formation 
rate shows that early peak in the case of HmNoRad (and 
HmNoDiff), but not in the cases which include radiation 
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Figure 3. The [Fe/H] metallicity distribution functions in the 
solar neighourhoods of the five simulations employed here. The 
bottom-right panel shows the MDF of the Milky Way's solar 
neigh bourhood, based upon two sub-samples of stars selected 
from iHolmberg et all J2009t) , as well as that for Local Group 
dwarf Fornax, from iKirby et alj ll201ll) (see text for details). In 
each panel, the overlaid curve is the best- fit single component 
Gaussian to the aforementioned MDF; the associated FWHM of 
said Gaussian is listed in the inset to each panel. 



energy - i.e., gas that forms stars (relatively) rapidly af- 
ter accretion does not mix as much, and hence the broader 
MDF. 

We next undertook a quantitative analysis of the MDFs 
shown in Fig [3] including a determination of the skewness, 
kurtosis, and widths at a range of inter-percentiles of the 
distributions. These determinations are listed in Table [2] As 
both skewness and kurtosis are highly sensitive to the pres- 
ence of outliers, we imposed a fairly standard 5<r clipping 
to the distributions. To mimic the observational uncertain- 
ties associated with the determination of individ ual stellar 
[Fe/H] abundances, after iFenner fc Gibson! l|2003l ). the 'the- 
oretical' MDFs shown in Fig [3] were convolved first with 
either a 0.1 dex Gaussian (to mimic the GCS uncertainties 
- Holmberg et al. 2009) or a 0.5 dex Gaussian (to mimic 
the uncertainties with the Fornax data - Kirby et al. 2011). 
In Table each column has two numbers; the first is the 
relevant metric, as measured on the MDF convolved with 
a 0.1 dex Gaussian, while the second (in brackets) is that 
measured on the MDF convolved with a 0.5 dex Gaussian. 
As the simulated MDFs are typically much broader than the 
GCS uncertainties, the impact of the 0.1 dex smoothing is 
minimal. 

As inferred from the above qualitative discussions of 
the MDF and the AMR (©, MDFs of the simulated solar 
neighbourhoods are all (save for HmNoMinShut, whose ex- 
ceedingly flat AMR results in the elimination of essentially 
all tails, positive or negative of the MDF's peak) more nega- 



tively skewed than that of the Milky Way's solar neighbour- 
hood (from both the volume-limited GCScut sample of stars, 
and the unrestricted GCS sample) and the sample from For- 
nax. It must be emphasised though that the typical 0.5 dex 
uncertainty associated with the determination of [Fe/H] for 
individual stars in Fornax means that broadening the simu- 
lated MDFs, with their typical dispersions of ~0.1 dex, by 
0.5 dex, 'washes out' much of our ability to compare and 
contrast the higher-order MDF metrics, and hence the anal- 
ysis which follows emphasises the differences between the 
simulated MDFs and that of the GCS. The 'tail' of stars 
to the negative side of the peak should not be associated 
immediately with the traditional 'G-dwarf problem', since 
these fully cosmological simulations relax the 'closed-box' 
framework which is the hallmark of this problem. Instead, 
as noted earlier, it is the tightly-correlated AMRs which are 
driving the large negative skewness values; these AMRs do 
not resemble that of the Milky Way's solar neighbourhood. 
The different SFHs are certainly part of the difference, but 
as noted earlier, both the fiducial llmChab and HmNoDiff 
show SFHs not dissimilar to the exponentially-declining one 
of the Milky Way, and the coordinated AMRs remain re- 
sponsible for the larger negative skewness in both cases. An 
analysis of the kurtosis values for each distribution are con- 
sistent with this picture. Specifically, the simulations' kur- 
tosis values are all larger than those of GCScut, and as noted 
in iJ21 large kurtosis values are driven in part by the pres- 
ence of a 'peaky' MDF, but more importantly, the impact 
of extended, 'heavy', tails. These tails (postive or negative) 
are driven by the coordinated AMRs and are reflected in the 
generally large values of kurtosis relative to the Milky Way's 
distribution. 

Alongside the skewness and kurtosis determinations, we 
present four measures of the shape of the MDF, through its 
dispersion, or width, at different amplitudes. This is done 
via the width of the inter-quartile range (IQR), inter-decile 
range (IDR), inter-centile range (ICR), and the inter-tenth- 
percentile range (ITPR)B 

The metrics associated with these width measures re- 
quire some comment in relation to the information provided 
by Fig [3] Specifically, the best-fit single Gaussian fits over- 
laid in each panel show that grossly speaking, the Milky 
Way's and Fornax's MDF are broader than those associ- 
ated with the simulations0 At first glance, the IQR, ITR, 
etc, measures listed in Table [2] appear counter to this result 
(which are all, essentially, larger than the values found for 
GCScut, for example). It is important to remember though 
that, much like the case for skewness and kurtosis, these 
measures of the breadth of the MDF are sensitive to the 
impact of outliers in the tails of the distribution. 

It is particularly useful to note the quantitative impact 
of the role of metal diffusion in setting the width of the 
MDF in tails of the distribution. For example, in the solar 
neighbourhood of the Milky Way, the range in metallicity 
between the bottom and top 0.1% of the stars is ~2 dex. 



17 The IQR corresponds to the difference in metallicity between 
the 25% lowest metallicity stars and the 25% higher metallicity 
stars; similarly, the IDR corresponds to the difference between 
the 10% lowest and 10% highest metallicity stars; etc. 

18 Save for HmNoDif f , as noted in [|3] 
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Simulation/Datasct 


Skewness 


Kurtosis 


IQR 


IDR 


ICR 


ITPR 


HmKroupa 


-1.84(- 


1.21) 


3.83(2.59) 


0.30(0.54) 


0.67(1.13) 


1.59(2.72) 


2.49(4.34) 


llmChab 


-1.56(- 


1.15) 


2.43(2.37) 


0.41(0.60) 


0.85(1.28) 


1.71(2.96) 


2.38(5.04) 


HmNoRad 


-1.13(- 


0.93) 


2.45(1.88) 


0.26(0.47) 


0.52(0.92) 


1.44(2.07) 


2.39(3.73) 


HmNoMinShut 


+0.47(- 


0.29) 


0.94(0.57) 


0.13(0.48) 


0.26(0.93) 


0.69(1.79) 


1.97(3.26) 


HmNoDiff 


-0.91(- 


1.29) 


0.91(2.32) 


0.96(1.25) 


1.85(2.44) 


3.49(5.18) 


5.06(8.03) 


GCS 


-0.61 




2.04 


0.23 


0.48 


1.26 


2.63 


GCScut 


-0.37 




0.78 


0.24 


0.45 


0.94 


1.43 


Fornax 


( 


-1.33) 


(3.58) 


(0.38) 


(2.25) 


(2.75) 


(2.85) 



Table 2. Prim ary MDF shape charact eristics for the solar neighbourhoods of the five simulations described here, the two sub-samples 
base d upon the [Holrnb erg et al l (|2009|) GCS empirical d ataset are as described in the text and data for the Fornax dwarf galaxy taken 
from lKirbv et al] [fcOlif T After iFenner fc Gibson! l|2003tl . the simulated MDFs were convolved with either a 0.1 dex Gaussian (left-most 
entry within each column) or a 0.5 dex Gaussian (right-most / bracketed entry within each column), to mimic the typical uncertainties 
associated with the [Fe/H] determinations in nature (the GCS in the case of the former, and Fornax in the case of the latter). Column 
(1): the name of the simulation or empirical dataset; Column (2): the skewness of the MDF (5<r clipping of outliers was imposed, to 
minimise their impact on the determination); Column (3): the kurtosis of the MDF, again with the adoption of 5c clipping; Columns 
(4) — (7): the interquartile (IQR), interdecile (IDR), intercentile (ICR), and inter-tenth-percentile (ITPR) for each MDF. 



For our simulation in which metal diffusion was neglected 
(llmNoDif f ), the corresponding width is ~5 dex - i.e., a 
factor of ~1000x greater than the other simulations with 
diffusion and that encountered in the Milk y Way, similar to 
what we found for other low diffusion runs (|Pilkington et all 
I2012D . 

After Casag rande et all (|201lh . we show in Fig [4] the 
MDF for the solar neighbourhood of one of our fiducial simu- 
lations ( HmKroupa), but now binned more finely in metallic- 
ity and colour-coded by age. Here, young stars correspond to 
those formed in the last 1 Gyr at redshift z—0; intermediate- 
age stars are those with ages between 5 and 7 Gyrs; old 
corresp onds to stars with ages g reater than 9 Gyrs. Using 
the GCS, ICasagrande et all (|2011 ) conclude that the younger 
stars have a narrower MDF that the older stars, consis- 
tent with our resu l ts (an d to be expected, given its AMR). 
ICasagrande et all l|201ll ) also found though that the loca- 
tions of the peaks associated with these old and young stars 
were at the same metallicity, which is not consistent with our 
simulations. Again, this is to be expected given the tightly- 
correlated AMRs of the simulations, relative to that of the 
Milky Way. 

While it may be the case that we are not capturing all 
of the relevant stellar migration physics within these sim- 
ulations (e.g., bars, spiral arms, resonances, etc.), there is 
radial migration occurring. That said, the radial gradients 
are shallow for these fiducial dwarfs (—0.01—0.02 dex/kpc, 
recalling Table and, as such, over the few kpcs of 'disc' 
associated with each simulated dwarf, systematic migration 
of metal-rich inner-disc stars outwards (and vice versa) has 
little impact on the position of the MDF 'sub-structure' (in 
which the young, intermediate, and old 'peaks' are offset by 



19 Flatter than the gradients seen in our work on the mas - 
sive galactic analogues to these dwarfs llPilkington et al.|[20l3) . 
consistent with the empirical work on gradients in dwarfs (e.g. 
ICarrera et "aT]|2008h . 
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Figure 4. The [Fe/H] MDF in the 'solar neighbourhood' of 
HmKroupa, split into three age intervals: young (black) defined 
as any star particle in the solar neighbourhood at redshift 2=0 
with an age less than 1 Gyr; intermediate (blue) defined as any 
star with an age between 5 and 7 Gyr; old (red) defined as any 
star with an age greater than 9 Gyr. 



~0.3— 0.5 dex from one another). Again, this is entirely con- 
sistent with the expected behaviour, based upon the AMR 
Fig 2. 

The central regions of our simulations show similar 
characteristics to those seen in the simulated solar neigh- 
bourhoods. Specifically, the [Fe/H] MDF of the 'bulge' (in- 
ner 2 kpc) shows a peak near [Fe/H]~— 0.5, with a number 
of sub-components at lower metallicity which correspond to 
progressively older and metal-poor populations (see Fig [5]). 
In spirit, such behaviour has been seen in t he M DF of the 
bulge of the Milky Way, where lfBensbv et alj l|201lf ) finds two 
populations, also separated comparably in age and metal- 
licity, to whic h they as s ociate seaprate formation scenar- 
ios. Similarly. iHill et al.l (2011) finds bulge sub-components 
within the MDF which they also separate into separate 
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Figure 5. The [Fe/H] MDF of the 'bulge' of llmKroupa; here, 
the bulge is simply defined as those stars located within 2 kpc of 
the galactic centre at 2=0. Alongside the full MDF (black line), 
sub-components based upon the age intervales noted in the inset 
are overdrawn. 



age, metallicity, and kinematic sub-structures, concluding 
the metal-poor component can be associated with an old 
spheroid, and the more metal-rich component can be asso- 
ciated with a longer timescale event (perhaps the evolution 
of the bar / psudeo-bulge). In our simulations, we see sys- 
tematic trends in age and kinematics for each metallicity 
sub-component of Fig 5, in the sense of the more metal-poor 
components being older and progressively less rotationally- 
supported, in exactly the manner one might predict from 
the AMR (©. It should be emphasised though that within 
the simulations, the behaviour of these age, metallicity, and 
kinematic 'sub-structure' in the bulge MDF is continuous, 
rather than showing any discrete transition from rotational 
support to anisotropic velocity support. 

Finally, we now examine in slightly more detail the be- 
haviour of the extreme metal-poor tails of the simulated 
MDFs (see Figs [6] and [TJl . In Fig [6] we show all stars be- 
yond the inner 3 kpc (and within 10 kpc), in order to min- 
imise the effect of the 'spheroid' stars in the analysis. We 
experimented, as before, with the impact of using a full 
kinematic decomposition between disc and spheroid stars, 
but again, for these dwarfs, the spatial cut alone is indis- 
tinguishable from the decomposed galaxy. In Fig [7] we only 
show those star particles lying within the previously defined 
'solar neighbourhoods' of each simulation. 

One additional curve is included in both figures (la- 
beled 109CH), that of the disc generated with the adap- 
tive mesh refinement code Ramses and described by 
ISanchez-Blazquez et al.l (|2009T ) . in which diffusion is han- 
dled 'naturally'. As noted previously, each o f the 11m series 
of simulations employ the lShen et alj (|2010T ) metal diffusion 
framework with a diffusion coefficient C=0.05, except for 
(obviously) HmNoDiff which assumes C=0.0. 

Each of the cumulative MDFs (Figs [5] and are nor- 
malised. In both cases, the normalisation occurs at the 
[Fe/H] corresponding to the metallicity of the lowest 1% 
of the stars (in terms of [Fe/H]). For plotting purposes, 
these are then aligned arbitrarily at [Fe/H]=+0.0, to show 
the relative distributions of extremely low- metallicity stars 



llmChab 
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Figure 6. The cumulative MDFs of the five 11m simu- 
lations: llmKroupa (black), HmNoRad (cyan), llmChab (blue), 
1 lmNoMinShut (gr een), and HmNoDiff (go l d), in addition to that 
of 109CH (orange: 'Sanchcz-Blazaucz ct al. ( 2009|) ) . For these six 
simulations, all stars lying within 3 and 10 kpc of their respective 
galactic centres are included in the analysis. The normalisation in 
each case is at the metallicity corresponding to that of the lowest 
1% (in terms of [Fe/H]) of the stars in each case. 



within each simulation and the empirical datasets. One 
could take a different approach and, say, normalise at (i) 
the same metallicity, (ii) the same amplitude, or (iii) the 
sa me number of st a rs. Fo r example, in our analysis of 
the iGovernato et al] (I2010T) b ulgeless dwarf galaxy simula- 
tions l|Pilkington et al.l 120121 ). we adopted (i), normalising 
all MDFs at [ Fe/H ]=-2.3. This was similar in spirit to 
ISchorck et ail (2009), who fixed the normalisations of the 
Milky Way halo and Local Group dwarf spheroidal MDFs 
to be unity at the metallicity corresponding to the lowest 
(in terms of [Fe/H]) ~100 stars in each. For distributions 
which peak at (potentially) very different metallicities, such 
normalisations can result in significant outliers which are 
not necessarily driven by any MDF 'tail'0 For our work 
here, while small quantitative differences exist depending 
upon the adopted normalisation, the qualitative results are 
robust regardless of the choice. 

What is immediately clear from even a cursory ex- 
amination of Fig [6] is that the relative distribution of ex- 
tremely metal-poor stars within all the simulations in which 
metal diffusion acts - i.e., all but HmNoDiff - are con- 
sistent with each other. This reflects graphically what we 
have commented upon earlier in relation to the tabulated 
ICR and ITPR values for the various MDFs (Table [2}. 
Specifically, the lack of metal diffusion within HmNoDiff 
drives its discrepant ICR and ITPR values (Tabled, and 
its outlier status in Fi g [6] When compared with Fig 5 of 
IPilkington et al.l |2012l ). one can see that the overly 'heavy' 
metal-poor tail to the MDF of HmNoDiff matches that en- 
countered in, for example, the low metal-diffusion simula- 



20 In the case of the analysis of lSchorck et al. (2009), the similar- 
ity of the positions of the peaks of the Milky Way halo and Local 
Group dSph MDFs meant that their analysis was robust against 
the choice of normalisation. 
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Figure 7. The cumulative MDFs of the analogous solar 
neighbourhoods associated the five 11m simulations: llmKroupa 
(black), HmNoRad (purple), llmChab (blue), llmNoMinShut 
(gree n), and llmNoDlf f (gold) , in ad dition to that of 109CH (or- 
ange: Sanchcz-Blazqucz ct al. (2009)). For these six simulations, 
the solar neighbourhood is defined spatially to include stars 
lying between 3 and 3.5 disc scalelengths from their respective 
galactic centres, and within 0.5 kpc of the mid-plane. The GCS 
and GCScut sub-samples described in ij2] are shown in red. The 
normalisation for each curve is as described for Fig [6] 

tions of lGovernato et alj (12010) Fl One fairly robust conclu- 
sion that can be drawn from Fig [6] is that the relative dis- 
tribution of extremely metal-poor stars is robust against the 
choice of feedback scheme; instead, diffusion plays a more 
important role in shaping this distribution. 

In some sense, the better 'statistics' afforded by Fig [6] 
provides a 'cleaner' picture than that seen when restrict- 
ing the analysis to just the 'solar neighbourhoods'[ff| For 
completeness though, in Fig [7] we also show the cumulative 
MDFs of the metal-poor tails for each dataset, normalised as 
in FigHJ] We should emphasise though that the small number 
of star particles in the 'bottom' 1% (in terms of metallicity) 
of the HmNoDiff, llmChab, and GCScut samples (~30 in 
each) make any interpretation susceptible to small-number 
statistics (and stochastic point-to-point 'fluctuations' which 
are 'averaged' over when employed the full disc, as in Fig [6]). 



5 SUMMARY 

Employing a suite of five simulations of an M33-scale late- 
type disc galaxy, each with the same assembly history, but 
with different prescriptions for stellar and supernovae feed- 
back, initial mass functions, metal diffusion, and supernova 
remnant cooling 'shut-off' period, we have analysed the re- 
sulting chemistry of the stellar populations, with a partic- 
ular focus on the metallicity distribution functions and the 
characteristics of the extreme metal-poor tail of said distri- 
butions. 

21 Demonstrating the quantitative power of the MDF to con- 
strain the magnitude of diffusion within SPH simulations of 
galaxy formation. 

22 And given the lack of any substantial gradient in the stellar 
populations for these dwarfs, the comparison is not invalid. 



In the context of the distribution of metals (in the sense 
of the higher-order moments of the resultings MDFs) within 
these discs, the impact of feedback and the IMF is more sub- 
tle than that of , for example, metal d iffusion. Employing a 
IChabrieJ (120031 ) IMF, rather than the lKroupa et al.1 (|l993T ) 
form adopted in our earlier work, does impact significantly 
on the resulting star formation history (and associated, re- 
duced, stellar mass fraction, resulting in remarakably close 
adherence to a wade range of empirical scaling relations - 
Brook et al. (2011c)). 

The star formation histories of the 'solar' neighbour- 
hoods associated with each simulation show exceedingly 
tight age-metallicity relations. In shape, these relations are 
akin to those pred icted by classical galactic chemical evo- 
lution models (e.g. IChiappini et al.ll200ll ; iFenner fc Gibson] 
2003), but bear somewhat less resemblance to that seen, 
for example, in the Milky Way's solar neighbourhood 
(Holmbcrg ct al. 200^). These correlated age-metallicity re- 
lations result inexorably in (negatively) skewed MDFs with 
large kurtosis values, when compared with the Milky Way. 
Star formation histories of dwarf irregulars, which quali- 
tatively speaking are a better match to those of the 11m 
series of simulations presented here, suggest though that 
somewhat steeper age-metallicity relat ions might eventuate 
in nature in these e nvironments (e.g. iDolphin et al.l 120031 : 
ISkillman et"ai] I2003T ). MDFs and AMRs of a com parable 
quality to that of the GCS (|Holmberg et al.ll2009h will be 
required to subtantively progress the field. 

An excess 'tail' of extremely metal-poor stars (amongst 
the bottom 0.1—1% of the most metal-poor stars) - 
~2— 3 dex below the peak of the MDF - exists in all of 
the simulations, as reflected in their inter-centile (ICR) and 
inter-tenth-of-a-percentile (ITPR) region measures. This tail 
is particularly problematic in simulations without metal dif- 
fusion (llmNoDiff) and those for which the diffusion co- 
efficient was set relativ ely low (e.g. iGovernato et aL I l2010l ; 
iPilkington et al.1 |2012| ). As demonstrated, the ICR and 
ITPR, in the absence of metal diffusion, can be ~30— 3000 x 
larger than that encountered in the Milky Way. 

We end with a re-statement of our initial caveat. 
The simulations presented here (particularly the fiducials, 
llmKroupa and llmChab) have been shown to be remark- 
ably consistent w ith a wide range of scaling relations 
l|Brook et al.ll2012l ). That said, their star formation histories 
are more akin to those of NGC 6822, Sextans A, WLM, and 
to some extent, the LMC (at least in the case of llmChab) 
- i.e., these systems are not 'clones' of the Milky Way. 
We have used the wonderful Geneva-Copenehagen Survey's 
wealth of data to generate empirical age-metallicity relations 
and metallicity distribution functions against which to com- 
pare, but exact one-to-one matches are not to be expected. 
That said, they do provide useful, hopefully generic, rela- 
tions against which to compare. In the future, we hope to 
extend our analysis to equally comprehensive datasets for 
the LMC, making use of, for exampl e, the data provid ed by 
the Vista Magellanic Cloud Survey (|Cioni et al.[|201ll ). 
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